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Ultrafast intersubband excitation of electrons in tunnell-coupled wells is studied depending on the 
structure parameters, the duration of the infrared pump and the detuning frequency. The temporal 
dependencies of the photoinduced concentration and dipole moment are obtained for two cases of 
transitions: from the single ground state to the tunnel-coupled excited states and from the tunnel- 
coupled states to the single excited state. The peculiarities of dephasing and population relaxation 
processes are also taken into account. The nonlinear regime of the response is also considered when 
the splitting energy between the tunnel-coupled levels is renormalized by the photoexcited electron 
concentration. The dependencies of the period and the amplitude of oscillations on the excitation 
pulse are presented with a description of the nonlinear oscillations damping. 



PACS numbers: 73.40.Gk, 78.47.+p 



-a 

C ■ I. INTRODUCTION 

o 

The coherent dynamics of electrons in heterostructures have been thoroughly examined during the past decade for 
the case of the interband ultrafast excitation by a near-infrared (IR) pulse (see Ref. [l| for review) . Recently, a mid- 
IR pump have been also employed for the treatment of the coherent dynamics of electrons under the intersubband 
excitation^. For example, a coherent transfer of electrons between tunnel-uncoupled states of a double quantum 
well (DQW) to the common excited state under mid-IR pump was considered in Ref. 0. Moreover, a new type of 
| semiconductor unipolar laser operating in the mid-infrared spectral region was demonstrated. This type of device 
("■-■«. . is based on a three-bound-state coupled DQW with a single-excited level and two coupled lower levels 4 . Thus, an 
investigation of the coherent dynamics in the tunnel- coupled DQWs under ultrafast mid-IR pump is now appropriate. 
In the present work we carry out the theory of the ultrafast response on the intersubband excitation between the 
tunnel-coupled states and the single level, which can be ground or excited. 

The study we will fulfill next is based on the quantum kinetic equation for the density matrix averaged over the pump 
frequency (see evaluation in Ref. 00). We will discuss the effects of the intersubband transition peculiarities by means 
of the intersubband generation rate. With this purpose we take into account the peculiarities of the intersubband 
. excitation for two cases: (A) when the electron transition occurs between the single-ground and the tunnel- coupled 
excited states, or (B) when the transition takes place from the tunnel-coupled states to the single-excited state. To 
illustrate these scenarios we have represented in Fig.l the band diagrams and the dispersion laws for two DQW samples 
O ■ of Ga As /AIq, 35 Gao. 65 As/GaAs, with the layer widths of 150/130/40 A and 150/20/120 A, corresponding to the cases 
(A) and (B), respectively. We have chosen the DQW structures in such a way that the energy separation between 
the coupled sub-levels, A T , is about 10 meV for both cases. In this context the population relaxation is controlled 
by the LO phonon emission^, while the dephasing of the tunnel-coupled states for the case (B) is determined by 
the quasi-elastic scattering. Since the interwell re-distribution of the charge appears under a relatively low pump 
intensity, we have considered both the second order response and the nonlinear regime of oscillations. Moreover, we 
will compare the present results with the corresponding ones to the interband excitation case. 

The paper is organized as follows. In Sec. II we derive the balance equations, which describe the coherent response 
of electrons in DQWs under the ultrafast intersubband excitation. In Sec. Ill we discuss the emerging quantum beats 
and the peculiarities of the coherent response under the finite duration excitation, stressing the differences between 
the cases of intersubband and interband excitation. Sec. IV contains the description of the nonlinear response. The 
conclusions and discussion of the approximations used are done in the last section. 
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II. BALANCE EQUATIONS 



The coherent dynamics of the electrons, when photoexcited by an ultra-short pulse, is described below in the 
framework of the second order response on the intersubband excitation. Performing the average over the peri od o f 
the radiation we obtain the quantum kinetic equation for the density matrix, p t , in the following form (see Refs. 1 51(1) : 



dp t 
dt 



[H,fit] = G t + i. 



(1) 



where H is the Hamiltonian of the DQWs under consideration, I sc is the collision integral, and Gt is the intersubband 
generation rate. When the electrons are excited by a transverse electric field E ± wt exp(— iojt) + c.c, with a frequency 
u) and a form- factor u>t, the generation rate is given by 
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Here A —> 0, the perturbation operator, Sht = (ie/uj)E ± v ± wt, is written through the transverse velocity operator v ± 
and p eq is the equilibrium density matrix when the second-order contributions to the response are taken into account. 

Neglecting the non-resonant mixing between the single and the tunnel-coupled levels we describe the system by the 
scalar distribution function, , where k = 0, ex correspond to the single electron state [ground |0) or excited \ex) 

state for the cases (A) or (B), respectively], and by the 2x2 matrix function f pt which describes the tunnel-coupled 
states \u) and \l) (upper and lower, respectively). Within the framework of the momentum representation, with the 
in-plane momentum p, Eq. (1) is transformed into: 

Pt Gg'+jW^tlP), 



dt 

gf f + r [hoQw, fpt] — G pt + I S c(ft\p), 



(3) 



where h DQW — (A/2)a z + Ta x is the matrix Hamiltonian of the tunnel- coupled states, A is the interlevel splitting 
energy, T is the tunnel matrix element, and a x . z are the Pauli matrices. Here the generation rates are different for the 
cases (A) and (B). Neglecting the overlap between |fc) and |Z) states, when (0|{i_i_|Z) — 0, and doing the straightforward 
calculations of Eq. (2), we obtain for the case (^4): 
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(4) 



where 9{e F — e p ) is the ground state equilibrium distribution for the zero temperature case, e F is the Fermi energy, 
and e p = p 2 /2m is the kinetic energy with the effective mass m. The dephasing time, T2, is introduced here instead 
of the A-parameter of Eq. (2) with the aim of describing a finite broadening of the intersubband transitions. For the 
case (B) we use (ex\v ± \l) ~ and the generation rate takes form: 
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(5) 



where p DQW is the equilibrium density matrix of the tunnel- coupled levels. The detuning frequency in Eqs. (4,5), 
Aw = u) — e /H, is evaluated through the energy difference between single and tunnel-coupled levels, e Q (see Fig.l). 
The remaining matrix elements in Eqs. (4,5) are calculated by using the matrix equalities: 
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Here Vl T — A T /h is the frequency of oscillations due to transitions between tunnel- coupled levels, A T = a/A 2 + (2T) 2 

and / e (±) = [6(e F -e- A/2) ± 6{e F -e + A/2)] /2. 

When doing the summation over the 2D momenta we introduce the population of the single level, N t — 
(2/ L 2 )^2 p F pt , and the 2x2 matrix of concentration (2/L 2 )^ p /p* — n t + (n* • &t), which is written through 
the scalar and vector components of the concentration, n t and n t . Due to the particle conservation law, N t +n t = n 2D 
with the total 2D concentration n 2D , the system (3) is transformed into the balance equations: 



^ - [L x n t ] + E(t) = G(t), 



(7) 



where S(t) = n"/ri for the case (A) or 5(f) = n t /ri for the case (B) and n" = n t + n\ . The vector £(i) is defined as 
£(£) = (0,0, n"/n) [case (A)] or £(£) = i>n t [case (£>)]. Here ti stands for the population relaxation time between 
single level and tunnel-coupled states, while the vector L = (2T/H, 0, A/ft) describes the dynamic properties of the 
tunnel-coupled electronic states. The relaxation matrix in the case (B), i>, is determined by the non-zero components 
(£)xz = {v)yy — Tq 1 , where the dephasing relaxation time, to, was introduced in Ref. 9 for the case of elastic 
scattering in DQWs. The generation rates G(t) and G(t) = [G x (t),G y (t),G(t)] are obtained from Eqs. (4-6) in the 
form: 
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(9) 



The photoinduced concentration in Eqs. (8,9) is determined as: 

7m 2D ( eE ± v ± 



(10) 



with the characteristic pulse duration t p and the characteristic velocities v 2 equal to |(0|wj_|/)| 2 or Kexlw^l?)! 2 for the 
cases (A) or (B), respectively. The coefficients a± in Eq. (8) are given by: a± = (1 ± An/n 2D )/2 [moreover An = 
for the DQW (A)] while, in Eq. (9), &± = 1 T A/A T for the case (A) and b± = (1 ± A/A T )(1 ± An/n 2D )/2 for the 
DQW (B), where An = p 2D A T . 

Next, taking into account the Coulomb renormalization of the tunnel-coupled levels, we have to replace h DQW in 
the matrix equation (3) by the Hartree-Fock Hamiltonian, h DQWl written in the form (see Refs. 10 and 11): 



= h L 
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(11) 



Here Q is the 3D wave vector, v Q is the Coulomb matrix element, and n Qt = Tr(p t e* Q r ) is the Fourier transform 
of the electron density. Further transformations lead to the balance equation (7) with the renormalized vector L t 
written through the level splitting energy 



A(t) 



A±^Z{nt 



(12) 



where Z is the distance between the centers of I- and r-QWs and e is the dielectric permittivity supposed to be 
uniform across the DQWs. The signs + and — in Eq. (12) correspond to the cases (A) and (B), respectively. The 
evaluation of A(t) coincides with that done for the DQW (A) in Rcf. 10. 



III. QUANTUM BEATS 

In this section we present a solution of the linear system of balance equations (7), neglecting the second addendum 
in Eq. (12), for the cases of short and finite pulse duration. Respecting the short-pulse approximation, if the pulse 
duration t p <C [Awl -1 ,!}" 1 , the generation rates [Eqs. (8) and (9)] take the forms: G x {t) ~ —a + (2T/A T )AfS p (t), 
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G y {t) ~ 0, and G(t) ~ b+j\fS p (t) with the S-\ike function: S p (t) = (2wt/ir) f®^ d,TWt+ T /T p . Thus, the photoinduced 
redistribution of the concentration can be written as the step-like function: n t = b + j\f f_ dt'S p (t') which is propor- 
tional to the step function 9(t) if t p — > 0. Since the photoinduced dipole moment is expressed through n\ , we obtain 
the z-component of n t in the form: 



n\ = 6(t)j\fUx>8 



cos (fl T t)}. (13) 



For the short-pulse approximation, the differences between the above-presented results and those corresponding 
to the case of the interband excitation (as considered in Ref. 5) are mainly attributable to the different character- 
istic concentrations and to the strong damping caused by the interband relaxation. Comparing Eq. (10) with the 
characteristic concentration for the interband excitation, N* [given by the Eq. (18) in Ref. Q, we obtain 



j\f 4n 2D ( E 



N* P* 2D {H/t p ) \E* Vs 



(14) 



where the interband excitation is characterized by the Kane velocity V, the gap e g , the reduced density of states 
and the field strength E* . If E ± ~ E*, and the pulse is not too short (t p ^lps), the ratio (14) is about 16 [case 
(A)] and 26 [case (B)) for the GaAlAs-ba,sed structures with a total 2D-concentration n 2D ~ 1.4 x 10 11 cm~ 2 and the 
dimensions used in Fig.l. Thus, the intersubband excitation appears to be more effective than the interband one. 

The response seems to be more complicated for the finite pulse duration case due to the peculiarities of the relaxation 
processes. We have used below the Gaussian form-factor, w t = exp[— (i/r p ) 2 /2], a semiempirical value of the damping 
To = 35 ps— , a dephasing time caused by the finite broadening of the intersubband transition T2 = 1 psA&i^ and an 
interband relaxation time due to LO phonons t\ = 3.5 ps^iI4. We consider first the evolution of the concentration. 
Fig. 2 shows the evolution of n t with the increase of the pulse duration T p f2 T /27r, for three detuning frequencies Auj = 0, 
Auj = Q T /2, and Auj — Q T [Figs. 2(a-c), respectively] and for the DQW (A). For DQW (B) the only difference 
is that the amplitude of the concentration n t is half of the corresponding to the structure (A) because, initially, 
there are two occupied levels in DQW (B). Therefore, we will pass by its interpretation, restricting ourselves to the 
case (A). One can see a new non-monotonic behavior in contrast to the one of the interband excitation case 5 . For 
< T p £l T /2ir < 1, n t /N behaves like in the interband case (corresponding to the short pulse context) with some 
type of oscillations superimposed. For t p VL t /2it > 1, nt/N these oscillations are strongly amplified around t = for 
Auj = and Auj = fi TS when the excited sublevel(s) is(are) not syntonized, before decaying. It should be noted 
that the excitation pulse is centered at t = 0. The number of oscillations depends on the pulse duration t p , as Figs. 
2(a,c) display. It is important to note that these oscillations have a period 2ir/Vt T , twice the n z t quantum beats 
period because such oscillations are controlled by the term Auj + Q T /2 and strongly influence the initial stages of nf. 
An exception takes place when one of the levels is syntonized, e.g., Aw = ±fi T /2. Then, the concentration shows a 
monotonous behavior with a growth rate similar to that of the interband pump [Fig. 2(6)]. Also visible in Fig. (2) is 
the exponential damping of the photoexcited electrons caused by the dephasing time, r 2 . 

Figs. 3 and 4 illustrate temporal evolution of the dipole moment, which is proportional to n\ , for different regions 
of parameters, t p , n = A/A T , and Auj. Figs. 3(a), 4(a) stand for the sample (A) and Figs. 3(6), 4(6) for the 
sample (B), respectively. The main difference between the finite pulse excitation and the short pulse excitation is 
the existence of two different regimes in the former event. When Auj = and r\ — [upper panels of Figs. 3(a,6)] 
the finite duration pulse produces a transition from a regime in which the electron density is mainly located in a 
well to two-well oscillations. This transition occurs when the pulse is switched off. The dipole moment exhibits the 
biggest oscillation amplitude while the pulse holds, then decaying due to relaxation until reaching the equilibrium 
after switching off the pulse. The balance situation is different for the two samples studied. In the first one the 
electronic redistribution between both wells are quickly reached, because the photoexcited electrons of the coupled 
levels decay to the ground state by means of the LO phonon emission. We must keep in mind that we are representing 
here the distribution n\ corresponding to the coupled excited levels. On the contrary, in the second sample, (B) , 
one can see the non-excited coupled levels. For this reason, the oscillations stay during some time until the electronic 
balance redistribution between the wells is reached because of the inter-subband dephasing relaxation. The time To 
for the last process is longer than that for the interband relaxation, t% (see numerical values above). Figs. 3 and 4 
show these features of the dipole moment in the cases of zero-phase shift (Q t t p — 4n) and 7r-phase shift (£I t t p = 57r) 
as indicated in figure captions. Fig. 3 has been calculated for r\ = 0, when the two tunnel-coupled states resonate 
and it corresponds to applied electric fields of 7 kV/cm (DQW (A)) and 2 kV/cm (DQW (B)), respectively. Fig. 
4 has been calculated for r\ = 0.7, out of the resonance of the tunnel-coupled levels. In this situation the electronic 
concentration mainly occupies the left well and the oscillation amplitude becomes quenched. 

The influence of the detuning frequency when 77 = can be explained as follows. If Auj — (upper panel of Fig. 3), 
a fast transfer of the electron density from the well in which electrons were initially created to the other well occurs. 
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For Auj = f2 T /2 (lower panel of Figs. 3), the electron density oscillates between coupled levels from the beginning of 
the excitation. Out of the resonance between the tunnel-coupled levels (77 ^ 0, Fig. 4) most of the electron density 
remains in the left well and the transfer doesn't become effective because of level decoupling. It is specially striking 
the practical disappearance of the oscillations when Auj = Q T /2. 

IV. NONLINEAR COHERENT RESPONSE 

Now we turn to the description of the nonlinear response. In order to do this we will take into account the Coulomb 
renormalization of the level splitting energy, when n t is governed by the nonlinear system of Eqs. (7), and L t is 
determined through Eq. (12). The characteristic concentration, TV, directly related to the pulse excitation density, is 
responsible for the nonlinearity. In order to get an effective Coulomb renormalization we have used TV > 2 x 10 10 cm~ 2 
(corresponding to an excitation energy density of about 10/1 J / cm 2 ) when the nonlinear response becomes noticeable. 

Figs. 5(a, b) show the evolution of the dipole moment, nf, corresponding to a characteristic concentration of 
TV* ~ 0.14n 2D (t p VL t /tt) , at the coupled-level resonance (77 = 0), zero-phase shift, and for structures (A) and (B), 
respectively. We should always keep in mind that TV depends on t 2 . Thus, for a fixed excitation energy, we have a 
different TV values for each pulse duration. The main result we can observe is that the oscillation period decreases 
and this is caused by a high TV value. This period also depends on the detuning frequency. As a consequence of this 
dependency, a slight Coulomb-induced dephasing appears between different t p and Auj cases. This behavior is more 
noticeable in the structure (B) than in (A) because of the relation jVire 2 Z/eT, which mainly determines Coulomb 
effects in Eq. (7) ( see Ref. |10|) . is greater in the former case for the same characteristic concentration because of the 
different values of v ± . Another feature induced by the Coulomb interaction occurs while the excitation pulse is acting 
on the samples. The term Auj + f2 T /2, which initially controls n t (and dipole moment oscillations), loses part of its 
importance and the masking of the intersubband oscillations diminishes. 

By comparing Fig. 5 with Fig. 3 one can see a slight displacement of the electronic concentration to the left QW 
caused by the above mentioned Coulomb renormalization when Aw = (upper panels). Once again the detuning 
frequency plays the main role in the oscillatory behavior, leading to a concentration, which is located in the left well, 
one order of magnitude higher for Aw = Q T /2 than for Auj — 0. Such a bearing is common for both samples studied. 

We have already shown (Fig. 4) that, being out of the resonance condition (e.g. 77 = 0.7), differences produced by 
the detuning frequency are small and this kind of behavior remains when the Coulomb renormalization is introduced 
[Figs. 6(a,6)]. However, there is a clear dissimilarity between structures (A) and (B). In the first sample the 
electronic concentration oscillates between the two wells from when the excitation pulse is switched on [Fig. 6(a)]. 
Such behavior is caused by a new situation of resonance at 7/ 7^ 0. To understand this point we must underline that 
the 77-values corresponding to resonance and off-resonance are strictly defined for the linear response. When the level 
renormalization is included resonance conditions vary and, hence, the electric fields to get them will also vary. In the 
other case, and for the same reason, electrons always prefer to stay mainly in the left QW [Fig. 6(b)]. These different 
behaviors are caused by the opposite sign in the expression for the Coulomb level splitting renormalization (Eq. 12). 
Finally, one can observe as a general bearing that the dipole moment oscillations are weak in the structure (^4). 
Furthermore, for both structures, the temporal evolution of the dipole moment loses its oscillatory behavior almost 
completely when Auj = Q T /2, the evolution depending essentially on the total concentration of excited electrons. 

V. CONCLUDING REMARKS 

Summarizing, we have described the coherent dynamics of electrons in DQWs taking into account the peculiarities 
of the intersubband excitation and relaxation for transitions between single and tunnel-coupled states. The temporal 
dependencies of the photoinduced concentration and the dipole moment are obtained both for the second order 
response and the nonlinear regime, when the splitting energy is renormalized by the photoexcited charge. 

Furthermore, we discuss the assumptions made. Both the tight-binding approximation for the description of the 
tunnel-coupled states and the use of the parabolic dispersion laws are valid for the DQWs under consideration. The 
simple relaxation time approach is also widely used for the description of similar structures. Applying the single- 
particle description of the high-frequency response we have neglected the Coulomb renormalization of the intersubband 
transitions due to depolarization and exchange effects, so that the nonlinear regime of the response under a not very 
low pump intensity may take place if A is not very big. On the other hand, we do not consider here the high-intensity 
pump case restricting ourselves to the inequality TV < n 2D when there is no Rabi oscillations. All these conditions are 
satisfied for the concentrations and intensities used in Sects. Ill and IV. 

To conclude, the peculiarities of coherent dynamics under the intersubband transitions of electrons described in 
sections III and IV are interesting in order to select effective conditions both for the THz emission, observed only under 



FIG. 1: Band diagrams and dispersion laws for the intersubband excitation of tunnel-coupled wells with single-ground (A) or 
single-excited (B) states. 



the interband excitation, and for the photoinduced concentration redistribution (see recent mid-IR measurements in 
a single QW**). It would also be interesting to verify scattering mechanisms by the use of this approach and to study 
the high-intensity pump, when an interplay between the nonlinear dynamics and Rabi oscillations appears. This case 
requires a special consideration. 
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FIG. 3: Temporal evolution of rit/Af for r\ — 0.. Figures 3a and 36 correspond to structures (A) and (B), respectively. Solid 
and dashed curves are plotted for t p — 1.76ps (zero-phase shift) and for t p — 2.2ps (7r-phase shift). Upper and lower panels 
correspond to Alo = and Alo = Q T /2. 
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FIG. 5: Nonlinear regime of the dipole moment nf/jV for r) = 0. Figures 5a and 56 correspond to DQWs (A) and (B), 
respectively. Solid and dashed curves are plotted for t p — 1.76ps (zero-phase shift) and for t p = 2.2ps (w-phase shift). Upper 
and lower panels correspond to Au — and Auj — f2 T /2. 




